Phantom crossing, equation-of-state singularities, 
and local gravity constraints in f{R) models 
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We identify the class of f{R) dark energy models which have a viable cosmology, i.e. a matter 
dominated epoch followed by a late-time acceleration. The deviation from a ACDM model {f = R — 
A) is quantified by the function m = Rf,RR/f,R- The matter epoch corresponds to m(r = —1) ~ +0 
(where r = —Rf,R/f) while the accelerated attractor exists in the region < m < 1. We find that 
the equation of state Wde of all such "viable" f{R) models exhibits two features: wde diverges 
at some redshift Zc and crosses the cosmological constant boundary ("phantom crossing") at a 
redshift zi, smaller than Zc- Using the observational data of Supernova la and Cosmic Microwave 
Background, we obtain the constraint m < 0(0.1) and we find that the phantom crossing could 
occur at Zf, > 1, i.e. within reach of observations. If we add local gravity constraints, the bound 
on m becomes very stringent, with m several orders of magnitude smaller than unity in the region 
whose density is much larger than the present cosmological density. The representative models that 
satisfy both cosmological and local gravity constraints take the asymptotic form ?n(r-) = C(— ) — 1)*" 
with p > 1 as r approaches —1. 
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I. INTRODUCTION 

Recent observations have continuously confirmed that 
about 70% of the present energy density of the universe 
consists of dark energy (DE) that leads to an acceler- 
ated expansion [H. The simplest DE scenario consistent 
with observations is the ACDM model in which DE is 
identified as a cosmological constant A. There arc many 
other attempts to explain the origin of DE, which can 
be broadly classified into two classes. The first class 
consists of modified gravity models in which gravity is 
modified from Einstein theory, whereas the second class 
corresponds to introducing a more or less exotic form of 
matter (such as a scalar field 0]) to explain the late-time 
acceleration. One of the main focus of current research is 
to find a departure from the ACDM model by confronting 
dynamical DE models with observations. 

In this paper we shall study cosmological and local 
gravity constraints on a class of modified gravity DE 
models 3 whose action is a general function f{R) in 
terms of a Ricci scalar R, i.e.. 



S = 



where 



ad 



(1) 



1/Afpj while G is a bare gravita- 



tional constant and Mpi is a reduced Planck mass (see 
also Refs. d, 0, Q). Here £m and £rad are the La- 
grangian densities of dust-like matter and radiation, re- 
spectively. Throughout this paper we shall focus on the 
metric-variational approach. See Refs. for the Palatini 
formalism of f{R) DE models. 

In Ref. Q it was shown that the models of the types 
f{R) = R- a/i?" in > 0) and f{R) = aR^ {n ^ 1) do 
not possess a standard matter epoch because of a large 
coupling between gravity and dark matter, even though 



a late-time acceleration can be realized (see also Refs. [1, 
[Tol . [HI). Extending the analysis to general f{R) cases, 
the paper [l^ has recently clarified the conditions under 
which f{R) DE models have a matter era followed by an 
accelerated expansion. However this does not necessarily 
mean that the models satisfying the conditions derived 
in [1^ can be consistent with observations. In this paper 
we constrain f{R) models that have a matter epoch prior 
to the acceleration from the observational data such as 
Supernova la (SNIa) and the sound horizon of Cosmic 
Microwave Background (CMB). 

The deviation from the ACDM model is quantified by 
a variable 



m = 



RR 



(2) 



where f^R = df/dR and f,RR = (Pf/dR"^. Note that 
the ACDM model corresponds to to — 0. As we will 
show in this paper, the quantity m can be constrained 
as TO < 0(0.1) throughout the matter and accelerated 
epochs from CMB and SNIa data. This limit still allows 
for interesting deviations from ACDM, in particular for 
the possibility to observe a phantom crossing and a sin- 
gularity at low redshifts in the equation of state of DE. 
However, as we will show, local gravity constraints (LGC) 
give a tight bound for to very much smaller than unity 
in a region whose density is much larger than the present 
cosmological density. Although this tight constraint ex- 
cludes most f{R) models (or renders them indistinguish- 
able from the ACDM model) , those proposed recently by 
Hu & Sawicki [isj and Starobinsky [l4| , which appeared 
after the initial submission of this article, are still viable 
and can be distinguished from ACDM. 
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II. VIABLE COSMOLOGICAL TRAJECTORIES 

111 the flat Friedmann-Robcrtson- Walker background 
with a scale factor a, the evolution equations in the 
metric- variational approach are given by 



3FH^ 



(Pm + Prad) + (F R - f)/2 - 3HF, (3) 



-2FH = K^[pn, + {4/3)p,,d]+F-HF, 



(4) 



where F = df/dR, H = a/ a, R ^ 6{2H^ + H), and a dot 
denotes a derivative in terms of cosmic time t. Note that 
we study the dynamics in the positive F branch. The 
energy densities of a non-relativistic matter and radiation 
satisfy the equations 



Prad + 4i7pi.ad = , 



(5) 
(6) 



respectively. 

In order to confront the models with SNIa observa- 
tions, it is convenient to write the equations as follows 

SFqH^ = K^(pdE + Pm + Prad) , 



2FoH 



(7) 

1^"^ [Pm + (4/3)p,.ad + Pde + Pde] , (8) 



where 



ajVde = {l/2)iFR- f)-3HF + 3H\Fo~F),i9) 
n^PDE = F + 2HF- {1/2){FR- f) 



{2H + 'SH^)iFo - F) ^ 



(10) 



Here the subscript "0" represents present values at the 
redsliift z = 0. By defining pde and puE in the above 
way, these satisfy the usual conservation equation 



Pde + (pde + Pde) = . 



(11) 



Then the DE equation of state (EOS) parameter, wde = 
Pde/pde, is directly related to the one obtained from 
observations [lH . 

Introducing the following variables 



F 



Xi 



HF' 
we obtain [1 



X2 



I 



a;3 



R K'^Piad 
, X4 = 



6H 



XlM = —l — X'i — ?>X2+x'\ — XiX'i+Xi, (12) 



m 
X1X3 



- ^2(2x3 -4-Xi) , 



-2x3(2:3-2) 



m 

-2.T3X4 + X1X4 , 



where iV = Ina, Xi^N = dxi/dN and 

^ dlogf _ Rf.HR 
dlogi? In ' 
^ ^ dlog/ ^ _Rj\R ^ X3 

d log i? / X2 



(13) 

(14) 
(15) 

(16) 
(17) 



Deriving i? as a function of r from Eq. pTj) . one can 
express m as a function of r = X3 /x2 and close the above 
system. Notice that defining fim = Pm/^FH^ one has 



\ — Xl — X2 — X3 ~ X^ . 



The DE equation of state is given by 

Pde 1 2x3 - 1 + x^y 



W-DE 



Pde 



3 1 - 2/(1 - a;i - a;2 - 0:3) ' 



(18) 



(19) 



where y = F/Fq. We note that the effective equation of 
state of the system is given by WcS = —(2x3 — l)/3. 

The analysis of the phase space has been performed in 
great detail in [l^ . Here we summarize the main results. 
In the absence of radiation {x^ = 0) we have six fixed 
points for the above system. A matter epoch can be 
realized for m sa and r « — 1 on the critical point 



Pm ■■ ixi,X2,X3) 



1 + Am 1 + Am 



3m 



1 + m' 2(1 + m)2' 2(1 -hm) 



which satisfies Wcff — — j^™^ . If m w +0 and m'{r) = 
dm/dr >— latr«— 1 the matter era corresponds to 
a saddle with a damped oscillation, whereas if m < a 
prolonged matter period is not realized because the real 
part of the eigenvalues of the critical point diverges (see 
also [i3l)- Note that the radiation point also exists in the 
region with m fa 0, which corresponds to a saddle p^ . 
Hence a viable cosmological trajectory starts around the 
radiation point with to sa 0, which is followed by the 
matter point with to w 0. 

When the trajectory passes through a standard matter 
era Pm with m w -1-0, then one sees that fim sa 1. On the 
other hand, the future asymptotic value of is always 
zero if the acceleration occurs (see below) . Therefore the 
denominator of wde, which in the absence of radiation 
can be written as 1 — Fflj^/Fo, goes from « 1 in the 
future to 1 — F/Fq in the deep matter era; this shows 
that if F increases toward the past 1 — F/Fq crosses zero 
and becomes negative and consequently wde passes nec- 
essarily through a singularity. We will show in the next 
section that this is indeed what happens. 

There are two stable fixed points leading to a late-time 
acceleration: 

(i) Pa ■■ (xi, X2, 0:3) = (0,-1,2), (20) 

(ii) Pb ■■ {xi,X2,X3) 
'2(1 -to) l-4m (l-4TO)(l-hm) 

1 -f- 2to 'TO(l-f2TO)' to(1 -f 2to) 

The effective EOS is given by Wcs ~ —1 for Pa and 
Weft = ^■,~^'?7o'"n for Pb- The former exists on the fine 

3m(l+2m) -° 

r — —2 and is stable for < 7ti < 1. The latter exists on 
the line m{r) = —r — 1 as is the case for the point Pm- 
There are several ranges of to that lead to a stable accel- 
eration, but only for (v^— 1)/2 < jti < 1 and m'{r) < —1 
one can have a transition from the saddle matter era to 
the accelerated epoch (in this case Woff > — 1) (see discus- 
sions in p^). Therefore, we have only two qualitatively 
different viable cases: 
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• Models that link Pm with Pa (Class A), 

• Models that link Pm with Pb (Class B). 

See Fig. [1] for an illustration. 

The cosmological dynamics of f{R) models can be 
well miderstood by considering m(r) curves in the (r, m) 
plane. The ACDM model, f{R) ^ R - A, corresponds 
to m = 0, in which case the trajectory is a straight line 
from Pm: {r,m) = (-1,0) to Pa: (r, m) = (-2,0). We 
now introduce three f{R) toy models that represent the 
two qualitatively different classes of cosmologies A and 
B. As Class A we define a minimal generalization of the 
ACDM model given by (model Al) 

/(i?) = (i?" - A)^ (21) 

This is characterized by the straight line 

m{r) = [{l-c)/c]r + b-l. (22) 

The existence of a saddle matter epoch requires the con- 
dition c > 1 and be 1. When c = 1 the trajectory is 
given by TO = 6 — 1 and is parallel to the ACDM line 
(to = 0). We also introduce models of the type (model 
A2) 

/(i?) = R- aR'' , a > 0, < n < 1 , (23) 

which satisfy m = n{l+r)/r and also fall into the Class A 
[ll|,[l2|. As representative of Class B, i.e., a cosmological 
evolution from Pm to Pb, we select the models (model 
Bl) 

TO(r) = -C(r + l)(r2 + ar + 5), (24) 

where C > 0. We require the conditions to'(— 1) = 
-C{l-a + b) > -1 and77i'(-2) = C{3a-b-8) < -1 for 
the transition from the matter era to the stable accelera- 
tion. In Fig. [T] we plot four trajectories corresponding to 
the various cases we presented above. These toy models 
are meant only to exemplify the classes of viable cosmolo- 
gies; it is clear that there are an infinite number of m(r) 
curves connecting Pm to Pa.b- However, most of our 
arguments below will apply to any of these curves. 

III. OBSERVATIONAL CONSTRAINTS 

For the viability of f{R) models we require that they 
satisfy three observational constraints: (i) CMB, (ii) 
SNIa and (iii) LGC. The first one comes from the an- 
gular size of the sound horizon defined by 

c.(z)dz / r-- dz 

where c'^{z) = l/[3(l + 3pf,/4p').)] is the adiabatic baryon- 
photon sound speed and z^ec — 1089 is the rcdshift at 
the decoupling time. This quantity Qg is related to the 
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Figure 1: 4 trajectories in the (r, m) plane. Each trajectory 
corresponds to (i) ACDM, (ii) f{R) = (i?*" - A)% (iii) f{R) = 
R - aR" with a > 0,0 < n < 1, and (iv) m(r) = -C{r + 
l)(r-^ + ar + b). Here Pm, Pa and Pb are matter, de-Sitter 
and non-phantom accelerated points, respectively. 



position of CMB acoustic peaks and has been constrained 
as &s = 0.5946 ± 0.0021 deg from the WMAP 3year 
data p^ . Since the eigenvalues of the Jacobian matrix 
for perturbations about the matter point Pm are given 
by 3(1 -I- to' (r)) and — 3/4 ± a/— 1/to, the matter epoch 
{m Ri +0) lasts for a long time as the tangent to'(— 1) 
approaches —1. If to' (—1) is close to —1 it is generally 
difficult to satisfy the CMB constraint. 

We have integrated the autonomous equations to ob- 
tain the present matter density rim"* = 0.28. Then 
we solved back the equations toward the past to get 
Qs — 0.5946 ± 0.0021. This is a trial and error proce- 
dure that is done by changing initial conditions in the 
radiation era. If the present radiation density satisfies 
the condition 5.0 x lO"'"^ < nfj^ < 2.5 x 10"'' together 
with the sound horizon constraint, we conclude that the 
models pass the CMB test. Then the straight line model 
Al given in (HI]), for example, is constrained to c < 3 and 
TO < 0.282. For larger c the contribution of dark energy 
is significant even in the matter-dominated epoch, thus 
incompatible with the CMB constraint. This comes from 
the fact that as to deviates from the coupling between 
dark energy and dark matter becomes significant. 

In what follows we shall consider the SNIa constraint 
under the situation where the CMB constraint is satis- 
fied. The viable cosmological trajectories are restricted 
to be in the range m > and r < 0. Since we are con- 
sidering the case of positive F, this translates into the 
conditions RF^r > and R/ f > 0. If i? changes sign, 
both to and r change signs simultaneously. Then it is suf- 
ficient to consider positive R, which gives Fji > 0. This 
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Figure 2: Evolution of the DE equation of state tode for the 
model f{R) = {R^'" - A)'= with parameters c = 1.01, 1.1, 1.8. 
As c approaches 1, the critical value Zc gets larger. In the 
limit c 1 (ACDM model) we have Zc — * oo. 

translates into the eondition F < Q when R decreases in 
time, which means that the variable F increases toward 
the past. Hence the DE equation of state exhibits a di- 
vergence at some redshift Zc > because pde changes 
sign. Since pde is negative provided > 1/2 (in which 
case the radiation, matter and accelerated epochs are in- 
cluded), we have wde < for z < Zc and wde > 
for z > Zc together with the singularity wde ^ — oo as 
z —Zc and wde — *■ +oo as z ^ +Zc- ^ This shows 
an interesting feature of f{R) models: whenever they 
are cosmologically acceptable they exhibit a singularity 
in Wde at some epoch in the past and, as a consequence, 
a crossing of the "phantom boundary" wue = ^1- Note 
that the models with F^n < (such as /(i?) = R — a/R" , 
n > 0) do not exhibit such peculiar behavior, but they 
are not cosmologically viable . 

In Fig. [2] we plot the evolution of wde for the model 
PT|) with be ^ 1 for three different values of c. The 
ACDM model (c = 1) corresponds to Zc ^ oo. As c de- 
viates from 1 , the critical redshift Zc gets smaller together 
with the increase of the present value of wde departing 
from —1. The phantom crossing (wde = —1) is realized 
at the redshift Zf, smaller than Zc. It is worth pointing 
out that the phantom crossing occurs from the region 
Wde < — 1 to the region wde > —1, which is different 
from quintom models of DE [l^. In Table I we show 
the values zi,, Zc, wde(-z = 0) and m(z = 0) with several 
different choices of c. We find that zi, is generally close 



We note that a divergence of iode can occur for a Dvali- 
Gabadadze-Porrati braneworld model [TtII . 



c 


Zb 


Zc 


w-de{z — 0) 


m{z = 0) 


1.01 


1.20 


6.77 


-0.996 


0.008 


1.1 


1.09 


3.33 


-0.952 


0.076 


1.5 


1.05 


2.55 


-0.818 


0.222 


1.8 


1.12 


2.52 


-0.766 


0.256 


2.3 


1.24 


2.61 


-0.705 


0.276 



Table I: The values of zt, Zc, wde(z = 0) and m{z — 0) for the 
model f{R) = {R^^"^ — K)'^ . The present epoch corresponds 
to = 0.28. 



Model 


Constraints 


Zc 


/ = {R^i" - ky 


m < 0.276, c < 2.3 


2.61 


f = R- aR" (0 < 7i< 1) 


m < 0.252, n < 0.7 


2.87 


m = -C{r + l){r + 2.1) 


m < 0.151, C < 0.5 


2.95 


m = -C(r + l)(r^+r + l) 


m < 0.295, 1/3 < C < 0.45 


2.40 



Table II: The constraint on the parameter m for several f{R) 
models coming from the SNIa constraint woEiz = 0) < —0.7 
and the CMB. We also show the value Zc corresponding to 
the maximum allowed ?n. 



to unity (unless iude today is extremely close to — 1), 
which is within the observational range of SNIa. This in- 
teresting feature could be employed to discriminate f{R) 
modified gravity from other dark energy models. The di- 
vergence of Wde occurs at a redshift larger than z = 2 (if 
one also imposes the CMB constraints) , so this is outside 
of the current observational range of SNIa. 

We can also use the criterion u;de(-z = 0) < —0.7 
for the compatibility with the SNIa data [131 ■ Then 
in the model Al we obtain the constraint c < 2.3 and 
m < 0.276, which is slightly stronger than the one ob- 
tained by the CMB. We have also carried out a numerical 
analysis for the other f{R) models A2 and Bl described 
in the previous section. In Fig. [3] the evolution of m is 
plotted as a function of z in the marginal cases satisfying 
both the CMB and SNIa requirements. The constraints 
on m for the models (Al) / = {R^/" - A)= ; (A2) / = 
R-aR"- (0 < n < 1), and (Bl) m = -C(r+ l)(r2 +r+ 1) 
are similar at all redshifts. Adopting instead a further 
model (A3) m = —C{r + l)(r -I- 2.1) (which also belongs 
to the class A) gives a tighter constraint near the present 
epoch. This is simply due to the fact that m decreases 
as r approaches —2 after having a maximum value of 
m at r = —1.55 in the (r, m) plane. The difference be- 
tween class A and B models is not significant provided 
the point Pb exists in the large m region close to 1 to 
ensure sufficient acceleration. In Table HIl we summarize 
the maximum values of m and the allowed model param- 
eters. As one can see, the parameter m is constrained 
to be m < 0.1-0.3 in all models we have considered. We 
have also checked that the slopes of the EOS, |dwDE/dz|, 
are smaller than «0.1 at present epoch and do not pro- 
vide better constraints than the ones obtained from the 
criterion wde{z = 0) < —0.7. 

As we have seen, the cosmological observations require 
TO to lie below 0.1-0.3 at all redshifts (up to the radiation 
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1 + z 

Figure 3: Evolution of the variable m in terms of z for the 
models: (Al) / = (i^^/'^-A)^ (A2) / = R~aR" {0 < n< 1), 
(A3) m = -C(r + l)(r + 2.1) and (Bl) m = -C{r + l){r'^ + 
r + 1). Each curve shows the maximal m that still satisfies 
SNIa and CMB constraints. 



epoch). For these values, the models display a phantom 
crossing and an equation of state singularity at a redshift 
of a few, making such f{R) theories quite intriguing from 
the observational point of view, especially since some SN 
analysis finds some evidence for a similar crossing (sec 

e.g. M)- ^ 

Let us now include the constraints from local gravity 
experiments. The Newtonian effective gravitational con- 
stant can be obtained under a weak-field approximation 
by considering a spherically symmetric body with a mass 
Mq , constant density p and a radius and a vanishing 
density (p = 0) outside the body. Using a linear per- 
turbation theory in the Minkowski background with a 
perturbation h^^^ and decomposing the function F into 
background and perturbation parts {F = Fq + SF), the 
effective gravitational potential at a distance £ from the 
center of the body is [ll [H, [H, [2^ 



1 



1 



where the mass AI is given by 



R 
3" 



f.B. 



B.R 



-Ml 



R 
3 m 



(1 - m) , 



(26) 



(27) 



If A/^ < the Yukawa correction e~^^^ in Eq. ([26|) is 
replaced by an oscillating function cos(|M|£), but this 
case is excluded experimentally. Hence the mass squared 
is required to be positive. 



As emphasized in Ref. [25| , the expression (pS)) is valid 
in the regime in which the linear approximation \5R\ <C 
Rq holds (here 5R is a perturbation in R with i?o being 



a background value). The condition for the validity of 
the linear approximation is given in Eq. (33) of Ref. [25| , 
which is equivalent to 



m(i?o) > $c 



(28) 



where m is defined in Eq. ^ and <I>c = GMq/tq is the 
gravitational potential at the surface of the body. Since 
the mass squared is estimated as ~ R/3m for 
m ^ 1, the condition ([28]) tends to be violated for large 
M. If we require Mi ^ 1 with i being the scale of the ex- 
periment (e.g. a scale of 1 mm for laboratory constraints 
or a scale of 1 AU for solar system constraints), we obtain 



(29) 



where Gn = G/Fq is the gravitational constant mea- 
sured at scales much larger than £. Note that we used 
the relation R w k'^p/Fq. This constraint is extremely 
stringent: assuming e.g. i = 1 AU and p = 10~^'^g/cm'^ 
(corresponding to the average density in the solar sys- 
tem) one gets m <C 10"^"^. Similar (or stronger) values 
are obtained for other experimental settings on the Earth 
or near other solar-system bodies. Combining Eqs. ([^5]) 
and we see that the condition for the applicability 
at a distance i from the center of a spherical structure of 
the linear perturbation approach is that 



£3p(£)> / p{i')l'^de' , 



(30) 



which, for most astrophysical bodies, is actually violated. 
In particular, since $c ~ 10~^ ^ 10~^ for the Earth or 
the Sun or other planetary bodies, we find that Eq. (1^5]) 
is actually grossly violated for local gravity experiments 
and we need to consider the non-linear regime. 

When the mass M is heavy, the system enters a non- 
linear stage in which a thin-shell develops inside the body 
through a chameleon mechanism [27| . In order to con- 
sider the chameleon effect in f{R) gravity, it is con- 
venient to transform to the Einstein frame by a con- 
formal transformation. Introducing a scalar field cf) as 
F = exp(-\/2/3K0), the potential in the Einstein frame is 
given by U = {RF — /)/2k^F^ with a constant coupling 
P = — 1/\/6 between the matter and the field (j) (sl. [2]|. 

In a spherically symmetric setting with an energy den- 
sity p, we obtain the following equation for the field (f> 



d^ 



2d(j) 
Id! 



dKff 



(31) 



where £ is the distance from the center of symmetry in 
the Einstein frame and 

Ueff (0) = V{cP) + ef'^^p* . (32) 

The energy density p* is defined by p* = e^^'^'^p, which 
is conserved in the Einstein frame 1271. 
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Let us consider a spherically symmetric body with an 
energy density p* = p*j^ inside the body (£ < Tq) and 
an energy density p* — Pb Pa outside the body 
{£ > r©). Then the effective potential ([32|) has two 
minima at (j) = <f)A and (j) = (pB satisfying the relations 
V4cj)A)+f3Ke''^'^^p\ = and +/3Ke'5''*«p|j = 0, 

respectively. The effective masses at the potential min- 



ima are defined by m 



2 — 
A 



and 



where the mass ttia is much heavier than the mass ms 
When a body has a thin shell, the solution to Eq. ((3T 
in the region £ > f© is approximately given by [l^, [27 



'47rMpi I 



where Mq = 47rr|p^/3 = 47rf|p^/3, 
PcS = 3/3— 



(33) 



(34) 



and $o = GMq/t-q. As long as the thin-shell con- 
dition ArQ/r-Q ^ 1 is satisfied, the effective coupling 
|/?off| becomes much smaller than unity. In this case 
the models can be consistent with the results of so- 
lar system experiments as well as equivalence principle 
experiments. For example, in the case of two identi- 
cal bodies with mass Mc, the potential energy associ- 
ated with the fifth force between the bodies is given 



by U{1) = 2(3i^{GM'Ji)e 



off 



< 



1221. 

10-3. 



The laboratory 
are satisfied for 



experiment constraints, 
Afo/fo < 1. 

In order to understand the condition under which the 
body has a thin shell, let us consider the model ((23)l with 
a ~ \Rl~"'. Note that Rc is not much different from the 
order of the present cosmological constant. In this case 
the Ricci scalar i?i at the de-Sitter point Pa satisfies the 
relation A = (i?i/i?c)^-"/(2 — n). In the region with a 
high density satisfying the relation Rc, the field (j>B 
and the parameter m{R) are approximately given by 



m{R) 



\/6. / Rc 

2 \K^PB 

An(l-n) 



M, 



pi , 



(35) 
(36) 



Using the fact that the Ricci scalar Rb in the region B is 
approximated by i?^ ~ k^Pb, we find that the thin-shell 
parameter is 



re 



1 



m{RB) 



2(1 -n) $0 



(37) 



This shows that, unless n is very close to 1, the thin-shell 
condition /^.tq/tq <^ 1 holds for 



(38) 



which is opposite to the condition p8l) for the valid- 
ity of the linear perturbation theory. When the non- 
linearity becomes important the body has a thin-shell. 
Since $0 ~ 10"^ and 10~^ for the Sun and the Earth 
respectively, the condition (j38p shows that the parame- 
ter m is very much smaller than unity in a high-density 
region where local gravity experiments are carried out 
{Rb > Rc). 

The current tightest constraint on the post-Newtonian 
parameter 7 in solar-system tests comes from Cassini 
tracking, which gives I7 - 1| < 2.3 x lO"'"^ [H, This 
translates into the bound I28l 



Are 



< 1.15 X 10"^ 



(39) 



Using this bound for Eq. ([37)1 with the value $© ~ 2.12 x 
10^^ of the Sun, we obtain 



Pi 



< 4.9 X 10"" 



(40) 



2-71 \pb , 

where pi = Ri/k?. Taking pi = 10^^^ g/cw? as the 
present cosmological density and pB = lO"^"' g/cvc? as 
dark and baryonic matter density in our galaxy, we ob- 
tain 

n<5xl0~*^. (41) 

Thus the model is very close to the ACDM model. The 
bound on n becomes even tighter if we take into account 
constraints from the equivalence principle [sj. 

Note that other /(i?) models we discussed in the pre- 
vious section also need to be very close to the ACDM 
model from the LGC. In such models the parameter m 
behaves as m = C(— r — 1) as r approaches —1 (here C 
is a positive constant). If we demand that the present 
value of m is of the order of 0.1 to find a deviation from 
the ACDM model, it is generally difficult to realize very 
small values of m around the region r w —1 to satisfy 
the constraint 

It is also worth mentioning that the model of Ref . [s^ , 
f{R) = R- XiRcexp{-R/X2Rc) with Ai, A2 > 0, which 
was explicitly introduced to satisfy the LGC without fine- 
tuned model parameters. When Ai « 1 the model passes 
the LGC for A2 < 10^. However, this model is not cos- 
mologically acceptable since it does not have a late-time 
accelerated attractor. We find in fact applying the cri- 
teria set forth in Ref. [l^ that (i) the de Sitter point 
Pa is not stable and (ii) the model does not have an 
intersection point with the line m ~ — r — 1 except for 
the point (r, m) = ( — 1,0), which implies that there are 
no additional accelerated attractors beside the unstable 
de-Sitter point. 

After the initial submission of this article, two viable 
models were independently proposed by (i) Hu & Sawicki 
[ist and by (ii) Starobinsky 



(i) f{R) = R-\R, 

(ii) f{R)^R 



{R/Rc 



l2ri 



1 



(i?/i?c)2" 

XRc [1- (l + R^/Riy 



(42) 
(43) 
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where n, A and Rc are positive constants. In such models 
the following relation holds in the region Rc [l^l : 

m(r) ~ C{-r - , (44) 

where C = 2n(2n+ 1)/A^". For larger n, m(r) decreases 
very rapidly as r approaches —1 to satisfy the LGC. 
Moreover, since m(r ~ —2) = C at the de-Sitter point 
Pa, it is possible to find a deviation from the ACDM 
model around the present epoch for C of the order of 
0.1. The detailed analysis about these models can be 
found in Refs. [l^, 0, l29l . [3^ . which showed that the 
models can be consistent with both cosmological and lo- 
cal gravity constraints for n > 2. Note that the model 
f{R) = R - APctanh (R/Rc) introduced in Ref. is 
also viable, which can be regarded as the limit n ~> oo in 
the models (j42|) and (j43|) . (see also Rcf. j3l| for a similar 
model) . 

IV. CONCLUSIONS 

We have shown that the variable m that characterizes 
the deviation from the ACDM model is constrained to 
be TO < 0(0.1) from the observational data of CMB and 
SNIa. We find that in general the f{R) models that are 
cosmologically acceptable exhibit a very peculiar behav- 
ior of the effective equation of state wde" this crosses in 
fact the phantom boundary and undergoes a singularity 
at a redshift of a few. If future observations will give 
the precise evolution of w'de hi the high-redshift range, 



it could be possible to detect the peculiar features of the 
f{R) models. This appears as an interesting way to dis- 
tinguish f{R) DE models from the ACDM cosmology. 

When we consider the local gravity constraints, wc find 
that the deviation parameter to is required to be very 
much smaller than unity in the high-density region where 
such experiments are carried out. In the f{R) models we 
studied in Sec. II, the parameter to has asymptotic be- 
haviour moc (— r — 1) asr approaches —1. If we demand 
that an appreciable deviation from the ACDM model can 
be found around the present epoch (to(z ~ 0) = 0(0.1)), 
we find that it is difficult to satisfy the LGC because 
TO does not decrease very rapidly in the region with a 
higher density. In such models viable cosmological tra- 
jectories satisfying all these constraints arc hardly distin- 
guishable from the ACDM model. However, the recently 
proposed models (j42|) and ([43|) can exhibit a deviation 
from the ACDM around the present epoch while satisfy- 
ing the LGC because the parameter to decreases rapidly 
in the higher-density region {R Rc). It would be of 
interest to place observational constraints on such mod- 
els including the data of matter power spectrum, SN la 
and weak lensing as well as LGC under the chameleon 
mechanism. 
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